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Abstract 

Mode coupling theory (MCT) has been successful in explaining the observed 
sequence of time relaxations in dense fluids. Previous expositions of this 
theory showing this sequence have required the existence of an ideal glass 
transition temperature Tq. Recent experiments show no evidence of Tq. We 
show here how the theory can be reformulated, in a fundamental way, such 
that one retains this sequence of relaxation behaviors but with a smooth tem- 
perature dependence and without any indication of Tq. The key ingredient in 
the reformulation is the inclusion of the metastable nature of the glass transi- 
tion problem through a coupling of the mass density to the defect density. A 
main result of our theory is that the exponents governing the sequence of time 
relaxations are weak functions of the temperature in contrast to the results 
from conventional MCT. 
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I. INTRODUCTION 



We show here how the condition of metastabihty in the glass transition problem can be 
used to establish the conditions necessary for a theoretical understanding of the experimen- 
tally observed elaborate sequence of time relaxation behaviors [|l| spread over many decades 
in time. This time sequence was originally predicted by mode coupling theory (MCT) [0,0 . 
We find a smooth temperature dependence for the relaxation without any evidence of an 
ideal glass transition temperature as hypothesized in conventional MCT. 

The mode coupling theory for the kinetics near the glass transition predicts a sequence 
of time behaviors which has received substantial verification through a steadily improving 
series of experiments in dense liquids This sequence, shown schematically in Fig.l for 

the density-density auto correlation function, begins, after those times associated with any 
fast microscopic processes, with a power-law regime / + Ait~"-. For times t c::^ Ta there is a 
cross over to the so-called von Schweidler regime where the correlation function decays as 
/ — A2t''. For times t > Tc one moves into the earliest stages of the primary or a relaxation 
which can be characterized by a stretched exponential behavior exp(— (t/r)^). There may 
then be a very long-time crossover to a final Debye or exponential decay regime. Within 
MCT the coefficients Ai, A2 and ^3 are parameters which depend on details of the model 
but the exponents a and b are related by a general nonlinear relationship |^ (discussed 
below) for which there is some experimental verification [^. For certain simplified models 
there is definite relationship between the exponent b and the stretching exponent (3. 
This relationship has been less well studied both theoretically and experimentally. The 
prediction of this sequence of time behaviors, which has substantial experimental support, 
must be considered a substantial achievement. 

Problems with MCT begin when one considers the temperature dependence gener- 
ated by the conventional assumption |^ that there is an ideal glass transition temperature 
To above and below which the dynamics are dramatically different. Indeed according to 
conventional MCT the von Schweidler and a-relaxation regimes are confined to the temper- 
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atures T > Tq. This sharp temperature dependence, as pointed out by Kim and Mazenko 
IP, is in disagreement with the high quahty experiments of Dixon et al. |T^. Furthermore, 
and again in disagreement with experiment, the conventional MCT finds that the sequence 
of relaxations must be confined to temperatures near Tq and that the associated exponents 
a, b and (3 should be temperature independent. The assumption in conventional MCT is 
similar to that associated with critical slowing down near a second order phase transition. 
Experimentally this picture does not hold. Dixon et al. ||10| find that the stretching ex- 
ponent (3 is weakly temperature dependent and that the qualitative relaxational behavior 
is the same for all temperatures spread around any reasonable choice for Tq. As pointed 
out by Kim and Mazenko, the data of Ref. 10 gives that (1 + 6)/(l + P) is temperature 
independent. Since P is temperature dependent the exponent b must also be temperature 
dependent. If the general MCT nonlinear relationship between a and b PJ21[] holds, then a 



must also be temperature dependent. MCT theory has even less to say about the remarkable 
scaling behavior found by Nagel and coworkers |T0[ . Kim and Mazenko p[ have shown that 
this scaling can be made compatible with MCT but at the expense of imposing additional 
ad hoc constraints on the theory. 

Our work here is to show how the MCT theory can be reformulated so as to be compat- 
ible with one's intuition that the substantial observed slowing down results from the fiuid 
system becoming metastable for temperatures below the melting temperature. The strong 
temperature dependence in the system goes into the r parameters which characterize the 



various frequency behaviors. It appears that these relaxation times show a power-law 
behavior as a function of temperature for higher temperatures but one expects such quan- 
tities to show an Arrenhius or Vogel-Fulcher temperature dependence for sufficiently low 
temperatures. 

Our analysis is based on a model that involves a number of simplifying assumptions. 
While most our final results are model dependent, our final conclusions may be quite general. 



The model we use is that due to Das and Mazenko , who developed fiuctuating nonlinear 



hydrodynamics for the mass density p coupled to the momentum density g, extended to 



include a coupling to an additional variable n which we associate with vacancy diffusion 
in the system. In the crystal we would be forced to include this variable, along with the 
Nambu-Goldstone modes corresponding to transverse phonons, in a rigorous treatment of 
the hydrodynamics of such systems ||13|. Here we assume that it is sensible to define such 
a variable in the dense but disordered state . One can think of this in terms of treating 



an order parameter in the disordered state. We will not need a microscopic definition of 
this variable. We assume that n is a diffusive variable with a diffusion coefficient F^. It is 
reasonable to assume that the time scale associated with n is very long, since one expects 
that defects must surmount an activation barrier in order to move. Thus we also assume 
that the driving effective free energy for n is a double well potential with the metastable 
defects associated with the higher energy well. We assume that the tunneling out of this 
well is facilitated by a coupling of n to the mass density in the effective free energy. We 
find self-consistently that the stretched dynamics can be associated with a defect potential 
with a very small activation barrier and a weak coupling between n and p. We show for 
the model we study that the vanishing of the barrier and the development of an inflexion 
point in the potential is a necessary condition for the slowing down observed experimentally. 
Apparently the system must build up the potential in which n moves. We also see that the 
notion of a 'below' the transition is not well defined from this point of view. 

II. THE MODEL 

The model we start with is the standard set of equations of fluctuating nonlinear hy- 
drodynamics for the set of variable p, g and n. We assume that p and n have the usual 
Poisson brackets of a scalar field with g and that the bare diffusion coefficients are just the 
bare shear and bulk viscosities rjo and (q for g and a bare diffusion coefficient Ty for n. The 
dynamics of this set is then driven by an effective free energy of the form: 

2p(x) 



d'^ PtT + Fu, (1) 



4 



where the leading term is the kinetic energy which follows from Galilean invariance. The 
term Fjj contains the dependence on p and n and will be discussed in detail below. The 
associated Langevin equations of motion are given then by 

I — . (.) 

^ = -/^^^^ - ^^^^ - E '^Apv^v,) - E u^v, + e. (3) 

| = -EV.(nV.) + r.V^^ + S, (4) 

where V = g/p and Lij{x) = ~r]o{^V iV j + 6ijV'^) — CoVjVj. The bare longitudinal viscosity 
is defined by Tq = Co + |^o- ©i and S in and are gaussian noises with variance 
depending on Lij and r„. To complete the specification of our model we must specify the 
potential Fu. Our simple choice is 

Fu = I d''x[j{5py + B{6p)n + h{n)], (5) 

where 6p = p — po and the defect potential is given by 

h{n) = eAiiir - w - + - i (6) 

where n is an average defect density in the metastable well in the absence of any coupling to p 
and e gives an over all scale for the potential and a measures the 'distance' from an inflexion 
point in the potential, or equivalently, the barrier height. In order to keep the analysis 
as simple as possible, we will assume in (j^) a very simple form for the density-dependent 
part of Fu which corresponds to a wavenumber independent structure factor. A closely 
related approximation is that the correlation function can be factored into a wavenumber 



dependent part and a frequency dependent part [|T^. This assumption works better than 
one might initially guess because, to a first approximation, the slowing down infiuences all 
wavenumbers. This is not a long wavelength or hydrodynamic approximation. 

The way in which one can include the energy variable in the problem has been discussed 
by Kim and Mazenko |TB|, but we neglect any such coupling here. The technical problems 



of treating the nonlinearities in this model using perturbative field theoretical methods flG 



has been discussed by Das and Mazenko [O] and clarified by Mazenko and Yeo |T^. The 



details of this analysis will be discussed elsewhere [0]. Here we briefly summarize the nature 
of the calculation and focus on the results. 

For the general set of equations described above with the assumption of the factoriza- 
tion of wavenumber and frequency dependent parts, one can show that the Laplace-Fourier 
transform of the density- density correlation function can be written in the general form 

^ / N J / N -2 + ^r(2;) 

where Qq is a microscopic 'phonon' frequency, 7(2;) is a long-time Das-Mazenko cutoff, and 
r(z) is the renormalized or physical longitudinal viscosity. Using the methods of Das and 



Mazenko |T2| one can calculate 7 and F in perturbation theory for almost any choice for 
the driving free energy. With this background and the assumption that current correlations 
decay much faster than p correlations. Das et al. indicated that the physical viscosity 
can be written in the general form 

r{z) = ro + nl dte'''H{<p{t)), (8) 

•/ 

where the mode coupling kernel can be written in the form 

N 

H{m = Y.(^nnt). (9) 

n=l 

The MCT analysis is carried out in terms of the coefficients c„. In the original 



Leutheussar model one had only the parameter C2- Das et al. [|20| showed how the 
development could be generalized to include an arbitrary set of parameters. We follow here 
the presentation of Kim and Mazenko A significant step forward was made by Gotze 
21| who realized that if one included a linear (ci) term in H then the von Schweidler and 
stretching behavior results. He included such a term on phenomenological grounds. Kim 
showed how such a term can be generated by coupling the mass density to an additional 



slower variable in the problem. One can also obtain the von Schweidler form and stretching, 
and effectively generate a ci term, if one includes the wavenumber in the analysis. This 
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has the physical interpretation that some band of wavenumbers is slower than others and 
stretches the faster modes. Some evidence for this picture was found in the simulations of 



Vails and Mazenko and the mode coupling calculations of Fuchs et al. ||2^ . Unfortunately 
an analytical treatment of these effects are not available. By including the wavenumber de- 
pendence in the development, as in the calculation of Das one can include such effects. 
Here we focus on the wavenumber independent model. 

The Das-Mazenko cutoff appearing in (|^) generates an exponential decay when it kicks 
in. We assume that the basic mechanism introduced by Sjogren and discussed by Kim 



and Mazenko [§ shifts the values of 7 to very small values. We assume for the rest of the 
discussion here, over the range of time scales we discuss, that this cutoff can be set to zero. 

If there was no coupling of p to n then our calculation is equivalent to that carried out 
by Das et al. pO[ and in more detail by Das and Mazenko [T^ and no linear term in H{t) 
is produced. The results in this case are equivalent to the Leutheusser model which 
has no stretching and no von Schweidler regime. The basic idea for introducing a Ci term 
in the analysis, as discussed by Kim and us elsewhere |TB[, is that the perturbation 



theory expansion for F will generate terms involving quadratic and higher order products 
of correlation functions. The correlation functions involving a current decay faster and 
those involving the vacancies slower than the density. Thus the vacancy-vacancy correlation 
function C„„ can be taken as a constant over the time range where the density-density 



correlation function Cpp is stretched [27]. Thus, for example, the quadratic form CnnCpp can 



be replaced by a constant times Cpp and generates an effective linear term in H. 

Metastability enters the problem since there is a significant time period over which the 
defects are trapped in a metastable state with some average density n. We restrict our 
perturbation analysis to this time regime. A key aspect of our analysis is that because n is 
coupled to p the stationary value of n depends on the local value of p determined by looking 
for the metastable minimum of F which is given by 

n* = n[l - — (^) - (1±^)(^)2(^)2 + . . .] (10) 
cry po a ay po 



where we have the dimensionless parameters x = Bp^n/Ap^ (a measure of the couphng 
between p and n) and y = en^ / Ap^. Our procedure is then to expand n about n* in the 
Langevin equations and then carry out perturbation theory in the nonhnear terms keeping 
terms up to = 2 in the mode couphng kerneL Self-consistently we find substantial slowing 
of the dynamics in the cases where x and a are very small. The smallness of a corresponds 
to a low barrier in h{n) while small x means that the potential is not greatly distorted by 
fiuctuations in p. The limit in which we obtain self-consistency is where x = Ca"^ with C of 
0(1). We then obtain, eliminating x in terms of C and a for small cr, that 

ci = i[-Ky + C) - + + SC)a] + 0{a') (11) 
C2 = e[l + % + 2(3y + 2C)a] + 0((t2), (12) 

where the formal expansion parameter ^ = kBTh!^ / (Qt^'^ Ap^) and A is a short-distance cut- 
off. The terms C3, C4, ■ ■ ■ are available if one goes to higher order in perturbation theory. 

We now demand that the parameters characterizing the defect well and the coupling are 
chosen to give the slow dynamics observed experimentally. This involves the mathematical 
machinery developed by MOT P,pT| but interpreted differently This analysis can be 
stated rather generally. If we define 

ao = (1 - f)VU) (13) 
= (1 - ffV'U) (14) 
\ = \{l-ffH"{f), (15) 

where V{f) = H{f) — //(I — /) and / has the physical interpretation as the metastable 
value of evaluated for the critical conditions ctq = en = which characterize the slow 
dynamics. In this case it can be shown that the decay exponents a and b are given by |21 



r(l-2a) r(l + 26) ^ ^ 

In the present case where we include only ci and C2 the critical conditions ai = can be 
used to eliminate the parameter / in equations (13-15) to obtain for small cxo, 

8 



Comparing (17) and (18) with (11) and (12) we obtain, assuming ctq — cr, four equations 
which can be solved to give C, y and A as functions of given by 

^ = i(l-w) (19) 

y = \{w^- 1) (20) 

and we have the imphcit equation for A as a function of ^ given by 

3(2A - 1) 

2[7 + 2A + V4A2 + 22A + 91]^' ^ ^ 

We can also express a in terms of ctq and ^: 

^ = "2[3y + 2C]eA(l-A)2''°- ^^^^ 
For I < A < 1, cr and ctq have the same sign. However, as can be seen from (6), only the 
absolute value |o"| matters in determining the critical condition, since the potential h{n) is 
'symmetric' |2^ around a = in the sense that both a > and cr < cases represent a 



double-well potential. 

III. DISCUSSION AND CONCLUSION 

Using the results of the previous section, we can express the exponents a and h and the 
amplitude / as weak functions of temperature, ^, as shown in Fig. 2. Note that we have been 
able, except for the temperature dependence of cr, to determine the form of the potential 
h{n) and its coupling to p as a function of temperature as represented by the parameter ^. 
At this stage it looks worthwhile to try and determine cr experimentally. 

Our analysis here has focused on the time regime prior to the c^-relaxation and we have 



not yet addressed the question of the scaling found by Dixon et al. |10]. Our detailed results 
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giving a, b and / as functions of ^ are model dependent and since our model is over-simplified 
should not be expected to apply to experiments in detail. However it is possible that our 
basic picture that a, b and / are weak functions of temperature and that the parameter a 
controls the slowing down may prove to be robust. 
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FIGURES 

FIG. 1. A schematic plot of the sequence of relaxation behaviors predicted by MCT; (a) 

power-law decay relaxation f+Ait^"", (b) von-Schweidler relaxation f—A2t^', (c) primary relaxation 
^36"^*/'^-*''; and (d) exponential relaxation e"'''*. 

FIG. 2. The exponent parameters a, b and the metastability parameter / as functions of the 
temperature represented by ^. The condition | < A < 1 restricts the range of ^ to 0.93 < ^ < 4. 
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